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Abstract. We present results for the 2m spectroscopic part of the Multi- 
Site Spectroscopic Telescope campaign, which took place in May/ June 2002. In 
order to perform an astcroseismological analysis on the multiperiodic pulsating 
subdwarf B star PG 1605+072 we used over 150 hours of time resolved spec- 
troscopy in order to search for and analyse line profile variations by using phase 
binning. We succeeded in finding variations in effective temperature and gravity 
for four modes. A pilot analysis using the BRUCE and KYLIE programs and 
assuming strong rotation and low inclination favours models with I = 1 or I = 2 
with m < 0. 



1. Introduction 

Subdwarf B stars are hot pre-white dwarfs (T e g = 22kK — 40kK) in the core- 
helium burning phase of their evolution (Heber, 1986). A new class of non-radial 
pulsators (EC 14026 or V361Hya stars) has been discovered (Kilkenny et al. 
1997). The pulsations are driven by an opacity bump caused by iron ionisation 
(Charpinet et al., 1996). PG1605 were discovered to be the largest amplitude 
pulsator amongst the V361Hya stars by Koen et al. (1998). More than 50 os- 
cillation periods have been found in its light curve (Kilkenny, 1999). Kawaler 
(1999) suggested a rapid rotation (v TOt = 130kms _1 ) to explain the frequency 
splitting. Indeed a high projected rotational velocity vsini = 39kms _1 has been 
measured (Heber et al., 1999). Time series spectroscopy revealed radial velocity 
variations (O'Toole et al., 2000) with more than 20 confirmed periods (O'Toole 
et al., 2005). In this paper we attempt to measure line profile variations and de- 
termine effective temperature and surface gravity variations in order to identify 
the individual pulsation modes. 



2. The MSST data and phasebinning 

In the 2m-part of the MSST campaign four observatories (Steward Observatory 
at Kitt Peak, ESO at La Silla, Siding Spring Observatory, NOT at La Palma) 
were involved. 9 nights of observation were followed by 12 more nights after a 
break of almost 3 weeks, resulting 10892 time resolved spectra (O'Toole et al., 
2005). We treat the data sets of each telescope separately. 
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In the first step we reduced spectra using the IRAF package. For each spectrum 
the Doppler shift, was removed. After that a continuum is fitted to the spectra. 
The line profiles are supposed to change if the star is pulsating, depending on 
its phase. The main reasons for those changes are of course temperature and 
density perturbations in the stellar atmosphere, besides the radial velocity vari- 
ations. To be able to detect these tiny variations, we determined the phase of 
the pulsation mode from the known periods for every single spectrum and coad- 
ded them accordingly. In this manner we got a complete pulsation cycle divided 
into twenty bins. Then we fitted LTE-model spectra to the binned spectra by 
performing a ^-minimization using the FITPROF program (Napiwotzki, 1999). 
So we were able to determine simultaneously the three atmospheric parameters 
(T e ff, log<7 and log He/H-r&tio) for every bin. Their variations for the strongest 
mode fl are shown in Fig. 1. The temperature amplitude is 1750K, while the 
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Figure 1 . The variations of the atmospheric values determined from the Steward 
data with error bars and sine fit are shown. In the upper panel temperature and 
gravity with sine fits are shown. The lower panel shows the He/H abundance and 
the temperature residuals. 

surface gravity variation is 0.16dex. The phase shift between the two fits is an 
indication for observing a nonradial mode. We also plotted the He/H abundance 
to make sure that it remains constant over a cycle. 

Fig. 2a shows temperature variations of the four strongest modes. Fig. 2a and 
Table 1 show the amplitude of the pulsations is decreasing with decreasing radial 
velocity amplitude of the mode. 

3. Modelling of line profile variations and mode identification 

In order to identify the modes we had to model various pulsation modes. We 
used the BRUCE and KYLIE routines by R.Townsend (1997). BRUCE calcu- 
lates an equilibrium surface grid of the star's envelope for a given set of stellar 
parameters. For any predefined set of quantum numbers it calculates effec- 
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Figure 2. a (left panel): the temperature variations for the four strongest modes 
determined from the Steward data with error bars and sine fit are shown. 
b (right panel): temperature and logg variations from synthetic spectra for the 
modes I = l,m = ±1 with sine fits and error bars are shown. 



Table 1. Amplitudes of temperature, surface gravity and radial velocity varia- 
tions for the four strongest modes. 
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tive temperatures and gravities for 40000 surface elements. KYLIE takes this 
grid and calculates synthetic spectra by interpolating in the same grid of model 
spectra used for the analysis of the observations. In the same way as with the 
observations, we used FITPROF to determine the atmospheric parameters out 
of the synthetic spectra from KYLIE for every phase bin. 

The required average atmospheric parameters (T e g = 32kK,logg = 5.25dex) 
have already been determined by Heber et al. (1999). Following Kawaler (1999) 
we choose a high rotational velocity of ISOkms" 1 and requiring a low inclination 
angle of 17° with a projected rotational velocity of vsini = 39kms _1 (Heber et 
al., 1999). 

In the last step we tried to compare the calculated amplitudes and phase shifts 
between T e g and log g with the measured ones. As there is almost no phase shift 
detected in the observations, modes with m < can be excluded. 
The analysis showed that the radial mode I = 0, m = has AT c g = 3355. 6K, 
which is much too high to be an observed mode, while the g-amplitude comes 
close to the observations. Further we discovered that for a fixed I always the 
mode with m = is the strongest. For all the modes with m < the tempera- 
ture and log g variations are almost in phase, while for the modes with m > 
there is a phase shift of typically ir. For the modes with I > 3 the variations 
became much too small to be detectable with this method (i.e. AT < 35K). The 
variations for I = l,m±l are shown in Figure 2b. While both modes have similar 
temperature amplitudes (m=+l: 539K, m=-l: 796K), the variations in surface 
gravity are drastically different: for m = +1 almost no variation (0.002dex) is 
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found, whereas for m = — 1 \ogg varies as much as 0.043dex. 
The T e ff amplitude predicted for the mode with I = 1, m = matches the ob- 
served fl pulsation (481.74s) quite well (s. Fig. 3). Although the surface gravity 
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Figure 3. Temperature trend for the strongest mode determined from the Stew- 
ard data and the adapted KYLIE model I = l,m = with error bars and sine fits 
are shown. 

amplitude is too small we regard this mode as the most likely. 
The second f2 (P = 475.61s) and third mode f3 (P = 475.76s) not only have 
almost the same period, but also the same amplitudes. In the modelling there 
were also only two modes I = l,m = — 1 and / = 2, m = — 1 which produced 
variations in that range. Therefore we believe to know the quantum numbers of 
f2 and f3 but cannot distinguish between them. 

Also for the fourth mode (P = 364.56s), there is actually only the mode with 

I = 2, m = —2 left, which produces variations in the observed range. 

The parameter range has to be further exploited to derive a consistent model. 
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